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Abstract 

We propose and investigate a spin-cell device which provides the necessary 

spin-motive force to drive a spin current for future spintronic circuits. Our 

spin-cell have four basic characteristics: (i) it has two poles so that a spin 

current flows in from one pole and out from the other pole, this way a complete 

spin-circuit can be established; (ii) it has a source of energy to drive the spin 

current; (iii) it maintains spin coherence so that a sizable spin current can 

be delivered; (iv) it drives a spin current without a charge current. The 

proposed spin-cell for spin current should be realizable using technologies 

presently available. 
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Traditional electronics is based on the flow of charge: the spin of electron is ignored. 
The emerging technology of spintronics will make the leap such that the flow of spin — in 
addition to charge, will be used for electronic applications [1,2]. A spin current is produced 
by the motion of spin-polarized electrons, therefore spin current is typically associated with 
spin-polarized charge current [1]. Nevertheless, if one can generate an ideal situation as 
shown in Fig. (la) where spin-up electrons move to the right while an equal number of 
spin-down electrons move to the left, then there will be no net charge current because 
7e = — where e/^, e/j are charge currents due to spin-up and spin-down electrons, 

respectively. There will be, however, a finite spin current: Is = |(-^t ~ h) where h is the 
reduced Planck constant. Considering the interesting and important future perspective of 
spin-current circuit, it is crucial to have a spin-cell that satisfies the four characteristics 
discussed in the Abstract and it produces the flow pattern of Fig. (la) [3]. In this paper we 
theoretically propose and analyze such a spin-cell. 

Our spin-cell is schematically shown in Fig. (lb). It consists of a double quantum-dot 
(QD) fabricated in two-dimensional electron gas (2DEG) with split gate technology, and 
each QD is contacted by an electrode. Note that no magnetic material is involved. The two 
QDs and their associated contacts to the electrodes serve as the "positive/negative" poles 
of the spin-cell. The two electrodes maintain the same electrochemical potential = I^'R 
{i.e. no bias voltage is applied on them). The size of the spin-cell structure is assumed to 
be within the spin coherence length which can be as long as many microns for 2DEG. We 
control the QD energy levels by gate voltages Vga where a = L,R indicates the left/right 
QD. Both QD levels are controlled by an overall gate voltage Vg, see gate arrangements 
in Fig. (lb). In order to distinguish spin-up electrons from spin-down electrons, a spatially 
non-uniform external magnetic field is applied to the two QDs-perpendicular to the QD 
plane. An extreme case of non-uniformity is Bji — —Bl, i.e. equal in value but opposite in 
direction. This particular magnetic field distribution is not necessary at all for the operation 
of our spin-cell, but it helps us to discuss its physics. Finally, the energy source of our spin- 
cell is provided by shining a microwave radiation with strengths Al/Ar for the left/right 
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QDs. Because, typically, the microwave frequency is far less than the plasma frequency of 
the material covering the QDs, the effect of the microwave field is to induce a high frequency 
potential variation AL/Rcoscut in the left/right QD and their leads. [4] When A^, ^ Ar, a 
time-dependent potential difference, Acosa;^ = (A^, — Ar) cos cut, exists between the two 
QDs. An a.c. electric field E(i) in the middle barrier is therefore estabhshed due to the 
microwave radiation (see Fig.(lc)). Then, electrons can absorb photons when they pass the 
middle barrier of the device. The establishment of E(i) across the two QDs is necessary for 
our spin-cell to work, here we use a non-uniform microwave radiation to achieve this effect 
as has already been carried out experimentally [5], but other possibilities also exist. 

Before we present theoretical and numerical results of the device in Fig. (lb), we first 
discuss why it works as a spin-cell. The physics is summarized in Fig.(lc). To be specific, 
let Br point to —z direction and to +z direction. Due to Zeeman effect, a spin-degenerate 
level cr on the right QD is now split into spin-down/up levels cr^ < cri- On the left QD, it is 

< ^Li- Electrons in the electrodes can now tunnel into the QD: on the right a spin-down 
electron is easier to tunnel because level e^^j^ is lower, while a spin-up electron is easier to 
tunnel into the left QD. Once levels eRi^e^ are occupied, the charging energies Ur,Ul of 
the two QDs push the other two levels e^^, eli to higher energies e^-f -|- Ur, eli + Ul, and 
the energy level positions indicated by the solid horizontal lines of Fig.(lc) are established. 
Next, the spin-down electron on the right QD can absorb a photon and make a transition 
to the level at e^^ -I- Ul on the left QD: afterwards it easily flows out to the left electrode 
because tii + Ul > A*l. This process is indicated as A—. Similarly the spin-up electron on 
the left QD flows out to the right electrode after absorption of a photon, indicated by A+. 
This way, driven by the potential variations of the QD induced by the microwave fleld, a 
spin-down electron flows to the left while a spin-up electron flows to the right of the spin-cell, 
and the continuation of the A—, A+ processes generates a DC spin current that flows from 
the left electrode, through the spin-cell, and out to the right electrode. Clearly, if the two 
processes are absolutely equivalent, there will be no charge current and only a spin current. 
Finally, since the spin-motive force is provided by a time-dependent change of the electronic 
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potential landscape of the QD, there is no spin flip mechanism and the spin current flowing 
through the spin-cell is conserved, i.e. Is,l = —Is,r — ^s- Our device then satisfies the four 
characteristics of a spin-cell discussed in the Abstract. 

The last paragraph discusses the operation principle of the spin-cell for spin current, but 
there are other interesting device details which can only be obtained by detailed theoretical 
and numerical analysis for which we now turn. The spin-cell of Fig. (lb) is described by the 
following Hamiltonian: [4,6] 

H= E + Wait) - {l/2)agnB^] dl^daa 

aa 

+ J2 Uadl^da^dl^dai + J2 i^^k + Wc,{t)] ali^^a^ka 

a kaa 
+ E [takaikadaa + H.c] + ^ [tcd{jRa + H.c] (l) 



kaa 



where dakai^ocka) and dl^^{daa) are creation (annihilation) operators in the electrode a and 
the dot a, respectively. The left and right QD includes a single energy level e^, but has 
spin index a and intradot Coulomb interaction Ua- To account for the magnetic field B, 
the left/right QD's single particle energy has a term —{l/2)agiJLBa, in which we have 
required a different magnetic field strength for the two QDs, i.e. B^ ^ Br. tc and 
Fq, = 2ti J2k \tak\'^S{^—^ak) describe the coupling strength between the two QDs, and between 
electrode a and its corresponding QD, respectively. The microwave irradiation is given by 
[4,6] Wa{t) — Aacosut and it produces an adiabatic change for the single particle energy. 
Here we permit the microwave field to irradiate the entire device including the electrodes, 
and we require a difference in the radiation strength Al ^ Ar. 

Our theoretical analysis of the spin-cell is based on standard Keldysh nonequilib- 
rium Green's function (NEGF) theory [4,6] which we briefiy outline here. First, we per- 
form an unitary transformation of the Hamiltonian with an unitary operator U{t) — 
exp \i Jq dt' ^ Wa{t')D >, where D = J2 CL^aka^aka + S d\,^daa- The Hamiltonian H is trans- 

l a J ka <y 

formed to the following form, 

-f^ = E ~ (^QfJ'Ba/'A dl^daa + E Uadl^d^^dl^dal 



kaa kaa 

La'- 



+E 



(2) 



where A = — A^. In (2), we take the last term which exphcitly depends on time t as 
the interacting part Hi, and the remaining part as Hg = H — Hj. The Green's function of 
Hoi g^{^), can be easily obtained with a decoupling approximation at the Hartree level [7], 

a"- (f) = ^gg + Ugl^aa /o^ 

yaaaK^) / _ n _l IP fe" + ?7 71 -V 
^aa)^aa ^ 2 °t\^aa ^ '-^a"'aaj 

where = e — Caa — Ua, eaa = — and Uaa is the time-averaged intradot 

electron occupation number at the state a in the a-QD which we solve self-consistently. 
It is worth mentioning that Qaaai^) Eq.(3) has two resonances: one is at Caa while its 
associated state at €aa is empty; the other resonance is at eaa + Ua while its associated 
state Eaa is occupied. Notice, in Hg the left part of the spin-cell {i.e. the left-lead and the 
left-QD) is not coupled with the right part of the spin-cell, therefore they are in equilibrium 
respectively. Hence the Keldysh Green's function Qaaai^) -^o can be solved from the 
fluctuation-dissipation theorem: g^aai^) — ~fa [9aaa{^) ~ 9aaa{^)] ■ With these preparations, 
the Green's function and of the total Hamiltonian H can be solved. In particular, 
we calculate Gl^^^it, t') = —i6(t — t') < {daaif), c^j^al^')} > by iterating the Dyson equation. 
In Fourier space, the Dyson equation can be reduced to [8,9] 

G^CT;mn(^) = ^a;mn{^) + ^^'a;mk{^)^'a;kn^^)^''r;nn{^) ^ 

k 

where G^,^^(e) = G^^„_^(e -|- ma;), and the quantity Gn{e) is the Fourier expansion of 
G{t,t'). [8] The retarded self-energy YT^.j^^{e) is the Fourier transform of S^(ii,t2) where 
^LWrihM) = ^^RLaitiM) = 5{h - h)tcexp[i dt' /^cosut'], and = E?j^ = 0. 

We obtain E£^^_^„(e) = '>y^La,nm{^) = tcJn-m{^)- The Green's function g;.„„(e) is: 
9ai3a;mn{^) = ^ap^mngaaai^ + "^^)- Then G^.^„(e) Can be solvcd from the above Dyson 
equation: [10] 



^aaa,mn{^) 



i^aaajmrn^ ■^c 



aaa,mn\ / o.aa,mrn\ ) ao.a,mniJ aaa,nn\ ) i 

where Aaammi^) = S \ tc\'^Jk-m \~) 9aaa-kki^) ■ Afterwards, the total Keldysh Green's func- 

k ^ 

tion G'^p^{t,t') = i < dp^{t')daa{t) > is easily obtained from the Keldysh equation. Finally, 
we obtain the time-averaged current in lead a from 

laa = < Iaa{t) >= -Im J {de / 2n)r ^{e) [G'<,,,oo(e) 

+ 2/a(6)G^«.,oo(e)] , (4) 

and the self-consistent equation for the intradot occupation number Uaa' — 
-tJ{de/27r)G<,,^,,{e). 

Fig. (2) shows the calculated charge current /g (in units of e) and the spin current Is 
(units of h/2) versus the gate voltage Vgn at different microwave frequency cu. 1^ shows a 
positive peak due to the A+ process and a negative peak by the A— process (see Fig.(lc)), 
but Is has two positive peaks. As we tune the gate voltage VgR, the right QD level is shifted 
so that when fiui — eli + Ul — e^j, the A— process occurs with high probability leading to 
a positive peak in Is and a negative peak in /g. On the other hand we get positive peaks in 
both 7e and Is when fiuj — er-^ + Ur — cli, for the the A+ process. 

The peak positions in Ig, Is due to the A-i: processes shift linearly with the microwave 
frequency a;, as shown by the dotted lines in Fig. (2). Eventually, at a special frequency 
indicated by A, i.e. when hcu* — er^ + Ur — —€^1 + 11^ — cr^, the two peaks overlap so 
that the net charge current Ig cancels exactly due to the cancellation of the A± processes, 
at the same time the spin current Is doubles its value. At this special frequency, the full 
operation of the spin-cell occurs so that a spin current is driven across the spin-cell, from 
left electrode to the right electrode, without a charge current. If we connect the spin-cell to 
complete an external circuit, a spin current will be driven and continue to flow across the 
spin-cell into the circuit. [11] On the other hand, if we let the two poles of the spin-cell open, 
although Is must be zero, a spin-motive force in the two poles of the spin cell will still be 
induced so that chemical potential iJ,ai 7^ jJ^ai- ^oi example, in the case of Fig.(lc), an open 
circuit will lead to /ili < /J'LI and /ir^ > /iri- 
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In the following we focus on the spin-cell operation by fixing gate voltage VgR — 0.45 
which is its value at point A of Fig. (2). We investigate le^h as functions of the overall 
gate potential Vg (Fig. (3a)), magnetic field q^Bl (Fig.(3b)), and frequency cu (Fig. (3c)). 
Different curves in Fig. (3) correspond to different microwave strength A = A^, — A^. In 
all situations /g ~ and we do not discuss it anymore. Fig. (3c) shows that Ig has several 
peaks and dips when we vary lu: the large peak indicated by A is the spin-cell operation 
discussed above, but peaks at C and D correspond to double- and triple-photon processes 
which connect the A± transitions of Fig.(lc). The dip at B originates from less probable 
transitions connecting levels indicated by the dashed lines of Fig.(lc), while the dip at E 
is its two-photon process. Now, fixing cu at a;*, i.e. at the spin-cell operation point A, the 
value of Is can be tuned by the overall gate voltage Vg as shown in Fig. (3a). However, Ig 
keeps large values for a wide range of Vg-. this range is in the Coulomb interaction scale 
U/e\ This is important, because in an experimental situation any background charge or 
environmental effect near the spin-cell may alter the overall potential, and Fig. (3a) shows 
that the spin-cell operation is not critically altered by this effect. When Vg becomes very 
large so that + Ul and + Ur is below the chemical potential or and is 
above /x, Ig diminishes because the A± processes can no longer occur (see Fig.(lc)). Finally, 
a very important result is shown in Fig. (3b), where we fixed q^Br = —0.4 while varying 
gjiBi at the spin-cell operation point [12] A. Fig. (3b) shows clearly that Ig increases with an 
increasing difference of {Bj^ — Br): Ig — identically when B^ — Br if Ul — Ur, or 7^ 
if Ul 7^ Ur. However, Fig. (3b) demonstrates that we only need a slight difference in Bl and 
Br, at a scale of the coupling constant F^, to generate a substantial Ig. The most important 
fact is that Bl and Br do not have to point to opposite directions which is experimentally 
difficult to do. In fact, if the two QD's are fabricated with different materials so that the 
gi— factors are different, one can actually use an uniform magnetic field throughout. 

The proposed spin-cell for spin current should be experimentally feasible using present 
technologies. First, the double-QD structures can and have been fabricated by several 
laboratories. Second, microwave assisted quantum transport measurements have recently 
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been reported [5,13,14]. In particular, the asymmetrical microwave radiation on double- 
QD device (i.e. ^ A^j) has already been carried out experimentally [5]. Third, the 
asymmetric magnetic field should be feasible as we have discussed above. If one takes / = 
ojjl'K — 50GHz, arranges the corresponding C/(~ 0.2meV, and fixes the temperature 

scale KbT and coupling r„ to be twenty times less than U as in typical QD experiments. 
I.e. ksT = lOOmK and T = lO/ieV, the corresponding magnetic field difference is {gii^Bi — 
Br) ~ r) \Bl — Br\ ~ 0.1 e/g' tesla. These QD parameters have already been realized by 
present technology. Finally, it is not difficult to show that by adjusting the gate voltages 
one can easily calibrate the spin-cell operating point [15]. 
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Appendix A 

In this appendix, we give a detailed derivation process of the unitary transformation in 
which the Hamiltonian (1) is transformed to Hamiltonian (2) in the manuscript. We take 
the unitary matrix U (t) as: 

U(t) = exp |z dt' E Wait') ai^^a^ka + E 4.^..) } (^1) 

Under this unitary transformation, the operator {Xa represents aaka and daa ) and 
transform into: 

U(t)XaUUt) = X^exp \-i f'dt'WJt')] (A2) 

L ^0 J 

U(t)XluUt) = Xlexp U f dt'WJt')] (A3) 

L Jo J 

Because unitary matrix U{t) is dependent of time the Hamiltonian H after the unitary 
transformation is: 

Hne^.{t) = U{t)H{t)U\t) - U{t)i^U\t) {AA) 



Using Eqs.(A2) and (A3), we easily obtain the new Hamiltonian (2) as given in the 
manuscript. We emphasize that this transform is done exactly and no approximation is 
involved. 

Appendix B 

In this Appendix, we present a detailed derivation for the retarded Green's function 
9aaai^) (Eq-(3) in the manuscript) of the Hamiltonian Hq. Here, we solve Qaaai^) ^y using 
the equation of motion method: 

e « A\B »'■=< {A,B} > + « [A,Ho]\B (El) 

where A and B represent arbitrary operators, and << A\B »^ is the Fourier expansion 
of the retarded Green's function —i9{t) < {A{t), B{0)} >. Because Hq is independent of 
time t, « A{t)\B{t') »^ only depends on the time difference t — t'. Appling equation of 
motion to the Green function g^^^^ —<< daa\dl^a have: 

1 4a +Ua « daad^daaldi^ »" {B2) 

k 

and 

Substitute Eq.(B3) into Eq.(B2), we obtain: 

(e - - S^) « daMl^ ! + d^^di^d^^\di„ »^ (54) 

ij. |2 

where EJ^ = '"+^0+ = "2^^^ is the self-energy function. 

At this point, if one takes a decoupling approximation for the new Green's function 
<< daadUdaa\dl^ >>'' m Eq.(B4), i.e., 

< < d^^di^daa 1 4a ^aa « C^aa 1 4a (^5) 

This is the mean field approximation. 
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Rather, we continue apply the equation of motion to << daadaadaa\daa 
have: 

(5 €cia << daa^d^^^da(f\d^^^ >> 'Tlaa 

+ 12^*ak« aakad^daaldi^ >>'^ 
k 

k 

+ Y.^lk«daadUa^ka\dL»' (BG) 

k 

In this equation, three new Green's functions << dakadaa-daaldaa- << 
daa^akadaaldl,^ >>^, and << daadlaO'aka\dl,^ »^ , emerge. We now take the decoupling 
approximation: 

<< aakadi^daa\di^ »''= Uaa « OaitaML >>'' 
<< daaalkadaa\dia >>''=<< daadl^ttakaldla {B7) 

Eq.(B6) reduces into: 

Then from Eqs.(B4) and (B8), we can solve the Green's function g^aai^) —« c^ao-Mao- »^ 
as: 

r / \ _ ^aa + UgJ^aa ( RQ\ 

This is the Eq.(3) in the manuscript. 

If one applies the equation of motion one more time, on the three new Green's functions 
of Eq.(B6), then takes a decoupling approximation, the Green function is obtained to be: 

9aaA^)- ,ip ,rr4 (yi-),yib)) ^^^^^ 

where the definition of Aaa-, and E^''^, are given in Phys. Rev. B 64,153306 (2001). 
This solution has a Kondo resonance and it has been applied in some previous work for 
investigating the Kondo effect. 
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Appendix C 

In this Appendix, we give a detailed derivation for the Green's function Ga;3o-,mn(^) 
the Hamiltonian H. First, from the Dyson equation, we have: 

^ LLa,mn 9LLa,mn ~^ ^ ^ ^ LR<7,mk^ RLa,kn9 LLa,nn (^■^) 

k 

k 

where we have suppressed the argument e. Substitute Eq.(C2) into Eq.(Cl), we obtain: 

^LLo-,mn ~ 9LLa,mn + 51 ^^LLa,mki^^LRa,kik2SRR(r,k2k2^^RLa,k2n9LLa,nn (^3) 

ki,k2 

Introducing notation An^^mni^), 

ARff,mni^) = ^^LRa,mk9RRa,kk^^RLa,kn (^4) 
k 

then the Eq.(C3) reduces to: 

^LLa,mn ~ 9LLa,mn ^ ^ LLa,mki^R'^Mn9LLa,nn 

ki 

~ 9LLa,mn 9LLa,mki^R^,kin9LLa,nn 
ki 

+ X] 9LLa,rnki^RcrMk29LLcT,k2k2^R<^Mn9LLa,nn + ••■ (^^) 

In the following we take an approximation as our previous work [Eq.(29) in Phys. Rev 
B 59,13126 (1999)]: 
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^R<T,kin9LLa,kiki9LLa,nn ~ ^RjMn^km 9LLa,nn i^^) 

This approximation is reasonable when hcu » 'max{ra,tc)- In the following we discuss 
the physical picture of this approximation. The Green's function G^io-mn ^ progapa- 
tor (of the Hamiltonian H) for an electron from the left dot to the left dot. The first 
term gLLa,mn describcs that the electron always stays in the left dot. The second term 

^ 9LLcTmki^Rt^Mn9LLann describcs the progapation process in which the electron from left 

ki 

dot transits to the right dot while absorbing ki — k photons, it then stays in the right dot 
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for some time, then goes back to the left dot emitting n — k photons. The higher order 
term describes multi-time transit processes between the left dot and the right dot. In our 
approximation {e.g. for the 2nd term), we take ki = n and neglect others. We emphasize 
that this term is indeed the most important. For example, if ki ^ n, we consider an initial 
electron in the left dot sitting at the resonance state {e.g. at e^). This electron tunnels to 
the right dot absorbing ki — k photons and goes back to the left dot emitting n — k photons. 
Then this electron has energy e^-f + {ki — n)hu!, which has a difference {ki — n)fiuj from the 
resonance state eL-f , and there is no state at energy ei] + {ki—n)huj. Therefore this processes 
{ki 7^ n) has very small probability to occur. In fact, when e is near a resonant state (e.g. 
^Li), Qaaai^) [Eq-(3) in the manuscript] can be rewritten as l/(e — + iS). Then the term 
in the left side of Eq. (C6) can be rewritten as: 



ARcr,kin 



(C7) 



(e — + kifiu + i5){e — + nUu + i8) 

where 5 is a small positive real number (~ F). Now it clearly shows that when ki ^ n, this 
term is always very small for any electron energy e [note that Aj^^^km has a small factor i^]. 
When ki — n, this term is large at e = eLi—nhcu. Therefore our approximation is reasonable 
if hcu » 'max{ra,tc). 

Under this approximation, we can easily solve G^^^^^^ from Eq.(C5) 

GLLcr,mn(^) = "7 T^T (^8) 

[9ll 

Then, G^ijo-^^n can be obtained immediately: 

<^W,mn(e) and G'-jiL {e) can also be solved in the same manner. 
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FIGURES 

FIG. 1. (a). Schematic diagram for a conductor which has a spin current with zero charge 
current; (b). Schematic diagram for the double quantum dot spin-cell; (c). Schematic plot for the 
spin-cell operation via photon-assisted tunneling processes indicated by A±. 

FIG. 2. The charge current /g and spin current Ig versus gate voltage VgR for different fre- 
quencies u. Different curves have been offset such that the vertical axis gives the frequency. Two 
dotted oblique lines A± indicate the position of the peaks. The parameters are: IJ-l = l^-R = 0, 
Tl = Tr = knT = 0.1, tc = 0.02, Ul = 1, Ur = 0.9, gf^B^ = 0.2, gf^BR = -0.4, VgL = 0.5, 
Vg = 0, and A/u = 1.0. 

FIG. 3. (a), (b), and (c) are /g and Ig versus gate voltage Vg, the magnetic field gfiBz, and 
frequency oj, respectively. u> = 1.25 in (a); 2u; = Ul + Ur + gfi{BL — Br) in (b). VgR = 0.45, 
and other parameters are the same as Fig.2. The solid, dotted, and dashed lines correspond to 
A/u> = 2.0, 1.0, and 0.5, respectively. Notice that the three curves of charge current overlap and 
they are all essentially zero. 
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